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ABSTRACT 



The gravitationally lensed Lyo' emitting galaxy, HCM 6A, detected by Hu et al. (2002) at z=6.56 behind the Abell 370 cluster was observed 
with the MAMBO-2 array of bolometers at 1.2 mm wavelength. The galaxy was not detected down to 1.08 mjy (3cr), but the depth of the 
observations and the lens amplification allow us to improve by approximately one order of magnitude previously published upper limits on far 
infrared emission of Lyo- emitting galaxies at this redshift. 

The following upper limits are derived from our observations assuming typical dust parameters: dust mass < 5.3 x 10^ Mg, IR luminosity 
< 2.1 X 10" Lo, and star formation rate < 35 Mq yr"' . The observed restframe UV-optical-IR spectral energy distribution (SED) of this galaxy 
is compatible with that of normal spiral galaxies or blue compact dwarf galaxies. SEDs of prototypical ULIRGs, such as Arp 220, are clearly 
excluded. Finally, we obtain an upper limit of S 2. 1 x 10"^ Mq yr"' Mpc"' for the dust-obscured SFR density of Lya selected galaxies at z ~ 
6.6. 
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^ 1. Introduction 



Nine of the 22 QSOs known to date at z~6 Jpaiil l2006t 
IWillott etal.''20O7) we re detected at mil l imeter wavelengths 
tPriddev et al,.2003: Bertoldi et al.ll2003at iRobson et alj|2004l: 
Wang etalj|2007h . The infrared luminosities implied by the 
millimeter fluxes are typically 10^^ L© and the dust masses are 
> 10^ M^. In addition, CO hnes were detected in J1148-I-5251 



dBertoldi et alJ l2003bt IWalter et al 



20031) revealing a large 
amount of dense molecular gas (2.2 x 10^" Mq) most proba- 
bly heated by a starburst rather than by an AGN. The Gunn- 
Peterson trough seen in the optical spectrum of these objects 
indicate that they are situated at the end of the epoch of reion- 
isation (Loeb & Bark ana 2001.). These observations imply that 
the process of enrichment of the ISM is relatively advanced 
in these objects and could have started at z S 8. This is in 
agreement with the results of W MAP which place the first on- 
set of star formation at z Sll jPage et al. 2006h . These ob- 
servations also put into question the processes responsible for 



the dust formation. Indeed, at the redshift of J1148H-5251 the 
Universe was only ~ 0.8 Gyr old and dust formation in qui- 
escent winds of low mass stars would not have been efficient 
enough to produce the observed masses. Possi ble sources of en- 
richment are SNII or pair instability SN (e.g. lTodini & Ferrara 
[200 1; Schneider et al. 2004), which can explain sever al obser- 
vations in the QSOs at z ~ 6 ( jVenkatesan et alj|2004l) . Quasar 
winds may also be a source of dust production. To better un- 
derstand the dust formation mechanisms and to obtain a more 
global view of star formation in the early Universe it is essen- 
tial to observe z > 6 non-quasar galaxies in the submillimeter 
range. To date submillimeter observations of only two ga laxie s 
at z > 6 have been reported in the literature (iWebb et al]l2()07l) . 
and each new observation brings important constraints to the 
models. 



The z = 6.56 galaxy HCM6A was discovered by lHu et al 
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(l2002h in a narrow band search for gravitationally lensed Lya 
emitters behind the galaxy cluster Abell 370 at ffjaooo - 
2'>39™54!73, cJjzooo = -1°33'32'.'3. Several follow-up obser- 
vations and analysis of this object have been published. In par- 
ticular, fr om an SED analys i s inclu ding constraints on its Lya 



emission. 



Schaerer & Pellol (120051) suggested that this galaxy 
could show a non-neghgible extinction, with A,, ~ 1. This re- 
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suit was later reinforced by the analysis of ICharv et alj (l2005b 
including also Spitzer observa tions, and more recently con- 
firmed bv lFinlator et alJ (l2007b . If true, the dust corrected star 
formation rate of HCM 6A is SFR ~ 1 1-41 Mq yr i.e. 2 to 
8 times higher than prev iously thought, and its luminosity is 
Lbo\ ~ (1 - 4) X 10" Lo ("Schae rer & Pell6ll2005l) in the range 
of luminous infrared galaxies (LIRGs). This high SFR implies 
the presence of large amounts of gas, and therefore dust. This 
makes this strongly lensed galaxy (magnification fj.=4.5 ac- 
cording to the model of Kneib et al. 1993) an interesting target 
to probe dust in very distant galaxies. 

Here we report new submillimetric observations with the 
MAMBO-2 bolometer at the IRAM 30m telescope. These 
observations were conducted in the frame of a larger pro- 
gram dedicated to near-IR and submillimeter observations of 
z > 6 g alaxies behind lensing clusters, as introduced by 



Schaerer et all (|2006|) . 

The paper is structured as follows: existing and new obser- 
vations of HCM 6A are described in Sect. 2. In Sect. 3 we es- 
timate the various physical parameters, and present the overall 
SEDs. Our main conclusions are summarised in Sect. 4. We as- 
sume a A-cosmology with Ho=70kms"' Mpc"', Qm=0.3 and 
Qa=0.7. 

2. Observations 

In Table [1] we have compiled the available observations for 
HCM 6A from the literature. All fluxes have been converted 
to micro Jansky; 3cr limits are given for non-detections. 



Table 1. Observations of HCM6A. 3cr limits are given for non- 
de tections. Note that a nearby source, CBK4, has been detected 
bv lCowie et al. I (l2002h at 850 fim. It is discussed in Sects. 12.1 
andim 



Filter/ 


wavelengtli 


tlux 


source 


instrument 


[jjm] 






V 


0.55 


<0.06 


Hu et al. (2002) 


R 


0.65 


0.02±0.01 




I 


0.86 


<0.06 




Z 


0.90 


0.15±0.11 




J 


1.25 


0.25+0.12 




H 


1.63 


0.40±0.15 




K 


2.12 


0.16+0.13 




IRAC 


3.6 


0.5+0.2 


Chary et al. (2006) 




4.5 


1.25+0.3 






5.8 


< 2.7 






8.0 


< 3. 




JCMT 


450 


< 30000. 


Small et al. (2002) 




850 


< 1700 


Cowie et al. (2002) 


MAMBO-2 


1200 


< 1080 


this paper 



2.1. Earlier observations 

The optical and near-I R data (up to the K' band) is taken from 
the discovery paper of lHu et alj (120021) . 



At longer wavelength the IRAC Spitzer fluxes from Chary 
et al. (2006) have been used. However, due to the presence of 
the ne arby source in the south-west already seen in images of 



Hu et al. (2002), the IRAC photometry may be somewhat con- 
taminated. For this reason and given the lower spatial resolu- 
tion, the MIPS 24 yum data cannot be used. 

The only references mentioning observations of Abe ll 370 
in the sub/milhm eter range are ICowie et al.l (l2002l) and 



Small et all (120021) . Both refer to SCUBA (JCMT) observations 
and the first reference only mentions the detection of a source 
(source No 4 in their Table 3, hereafter CBK4) with a flux of 
2. 17+0.57 mJy at 850yum and located ~ 10" to the West of 
HCM6A. This emission likely originates from a foreground 
galaxy of the cluster ~ 5" to the South- West of HCM6A. In 
addition this offset is larger than typical SCUBA astrometric 
uncertainties (~4"). However, the astrometry becomes less ac- 
curate for weaker sources and we note that for another weak 



source in the field the positions given by ICowie et al.l (l2002h 
and Small et al. (2002) could differ by 10". Indeed, the source 
#2 oflCowie etalJ(l2q02h detected at 5 cr might also be detected 
bv lSmail et alj ( |2002|) but offset by -10" (the source J02400- 
0134). Hence, it cannot be definitely excluded that the emis- 
sion of CBK4 comes from HC M 6A. The galaxy w ould be at 



the border of the field covered by' Smail et al. ( 200 2*) but its flux 



would be below their detection threshold. In short, except oth- 
erwise stated, we retain the upper limits from SCUBA in Table 
[T]for HCM 6A (for the 850 /urn upper limit we take 3 times the 
uncertainty quoted for the flux of CBK4, i.e. 1.7 mJy). 

2.2. New MAMBO-2 observations at 1.2 mm 

The observations were carried out at the IRAM 3 0m Telescope 
with the MAMBO-2 array of 117 bolometers dKrevsa et al. 
2002) during the pool session in 2004/2005. It has a half- 
power spectral bandwidth of 80 GHz centered on ~250GHz 
(1.2 mm). The beam size is 11", thus, any emission origi- 
nating from the foreground galaxy of the cluster ~ 5" to the 
South- West of HCM6A would be dimmed by a factor of ~2. 
The excellent weather conditions allowed us to reach an rms 
of 0.36 mJy in 4.2 hours of integration in ON/OFF mode. The 
galaxy was not detected, the flux measured is -0.03 +0.36 mJy. 
Our non-detection corresponds to a 3cr upper limit of 1 .08 mJy. 

3. Results 

The main derived quantities for HCM 6A are summarised in 
Table|2] We now discuss them one by one. 

3.1. Dust mass 

The upper limit on the MAMBO-2 flux allows us to compute 
an upper limit for the mass of dust using the fomula: 



(1) 



where 5 y is the flux upper limit, fi is the gravitational magni- 
fication factor, Dl is the luminosity distance, Ky^ is the mass 
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Table 2. Derived quantities for HCM 6A assuming different 
dust temperatures Tj, a magnification i-i-A-.S, a dust mass ab- 
sorption coefficient /('i25=1.875m^kg"', a dust emissivity in- 
dex y6=1.5, and a high frequency dust spectral index a-2.9. 
All limits are 3cr upper limits based on the 1.2 mm flux up- 
per limit for Td^lS and 36 K and the 850;um flux upper 
li mit for rd =54K. For comparison, quantities estimated by 
Schaerer & Pella (l2005h are given at the bottom of the table. 
Note, however, that the latter SFR estimates need to be multi- 
plied by 2.55 to account for different lower mass limits of the 
Salpeter IMF (Miow = 1 M© instead of 0.1 assumed for the 
SFR(IR) determination). 





Ti = 18 


Ti = 36 


7-d = 54 


Dust mass [Mq] 


< 7.0 X 10** 


< 5.3 X 10' 


< 1.5 X 10' 


i'FIR [Lq] 


< 6.4 X 10'" 


< 2.1 X 10" 


< 5.4 X 10" 


SFR(IR) [Moyr-'] 


< 11 


< 35 


< 87 
mil 



!^bol 

SFRuv(A,, 
SFRLy„(A,, 
SFRuv(A,, 
SFRLyff(A,, 



: l)[Moyr- 
= l)[Moyr- 
0) [Moyr- 



0) [Meyr-'] 



11-41 
7-12 
5.6 
0.4-0.8 



a. 



4.5 






\^ SCUDA©450/im 




4.0 


- T=54 


K 






3,5 






SCVBl 


MAMB0iai.2mm - 


3.0 


T=38 








2.5 










2.0 


T=18 


K 






1.5 











1.6 1.8 2.0 2.2 2.4 

Log(A„.,[^m]) 

Fig. 1. Far infrared SEDs for dust temperatures of 54 K 
(dashed), 36 K (soUd) and 18 K (dotted) with yS=1.5, a^2.9 and 
such that f i.2mm-l-l mJy. The triangle s hows the possible de- 
tection at 850 yum by Cowie et aU (l2002l) . the arrows show the 
3cr upper limits as summarized in Table[T] 



The dust mass estimate is highly sensitive to the tem- 
perature assumed, the upper limit becomes 2.0 x 10^ Mq for 
rd=54K and 7.0 x 10** Mq for rd=18 K. In addition, the dust 
at high redshift is thought to result mainly from condensa- 
tion of grains in SN ejecta and its emissivity might differ 
from that of the local Universe. With the emissivity introduced 
by Bianchi & Schneider (2007) for dust formed in SN ejecta 
we obtain 65% lower values for the upper limits on the dust 
mass. Uncertainties i n the em issivity index must also be con- 
sidered. According to lYang & Phillips (2007) the median value 
for LIRGs is y8=1.6+0.3. Taking the higher value, i.e. j6=1.9, 
increases the dust mass upper limits by 10%. However, despite 
the uncertainties, this upper limit definitely excludes the dust 
masses found in the high-z quasars and submillimeter galax- 
ies (SMGs). It is closer to the mass of dust observed in local 
Universe galaxies dVlahakis e t al. 2005). 



If the emission seen in the 850 fim map by ICowie et al 



(120021) is real and if the flux effectively originates from 
HCM6A (cf. above discussion about CBK4), it would imply 
a dust mass Md=6.5 x 10^ M©, assuming - 36 K. Taking 
into account the uncertainties this is compatible with our upper 
limit at the same temperature. Our upper limit would however 
imply that the dust cannot be colder than 36 K (see Fig.[T]). If, 
as is most likely, the emission seen in their map does not origi- 
nate from HCM6A, our 1 .2 mm dust mass upper limit improves 
their 850 /vm upper limit for dust colder than 36 K. For a higher 
dust temperature their upper limit is more stringent, e.g. at 54 K 
it would imply Md<1.5 x 10^ Mq. 

In the following and as summarized in Table|2]we allow for 
an uncertainty of 50% on the dust temperature by considering 
the three values, 18, 36 and 54 K. When refering to the 2 low- 
est temperatures we implicitly use our 1 .2 mm flux upper limit 
and f or the highest temp erature we use the 850 //m upper limit 
from Cowie et al.l (l2002h . The errors seem to be dominated by 
the uncertainty on the temperature r ather than the emissivit y 
index. In addition, if, as proposed bv lYang & PhillipsI (l2007h . 
the emissivity index is negatively coiTelated to the tempera- 
ture, then the variations of both parameters would compensate. 
To simplify the error analysis we therefore consider a range in 
temperature only and we assume that the resulting errors are 
good approximations to the errors due to all dust parameter un- 
certainties. 



3.2. Far infrared luminosity 



absorption coefficient of the dust at Vr, and By^iJS) is the in- 
tensity of a blackbody at Vr assuming isothermal emission 
from dust grains at tem perature Ta . To compare with the 



upper limits reported by IWebb et al.l (120071) we also assume 
''■i25//m=l-875 m^ kg"' and Ky oc \fi with the emissivity index 
jS=1.5 and we take rd=36 K. This leads to a 3cr upper limit of 
Md=5.3 X 10^ Mq. This upper limit is o ne order of magnitude 
lower than the two upper limits given by Webb et alJ (2007) for 
two z~6.5 Lya emitters. This significant improvement in the 
dust mass upper limit of a Lya emitter at z >6 results from the 
high sensitivity of the observations and the strong gravitational 
amplification of the source. 



As lWebb etall ( l2007h we compute the total far infrared (FIR) 
luminosity assuming a modified black-body spectral energy 
distribution (SED) with a dust temperature of 36 K and a 
dust emissivity of y6=1.5 and replacing the Wien's expo- 
nential decrease shortward of 53 /im by a power law with 
spect ral index a- 2.9 to better agree with observed SEDs 



(Blain et al.l 120031) . Our 1.2 mm 3cr flux upper limit corre- 



sponds to Lfir=1-9x10 Lq, which is one order of magnitude 
lower than the FIR luminosities usually measured in the SMGs 
co mparable to the local Univers e ultraluminous infrared galax- 
ies dChapman et alj|2004l2005l) . 
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This leads to an upper limit on the Lpjr/Luv ratio of ~2. 
This is clearly within the range o f values exhibited by Lyman - 
break galaxies of LpiR /L.Tv <10 dAdelberger & Steidell l2000l). 
Using the fits of Burg arella et alT l 2005b our Lftr/Litv limit 
translates into a dust attenuation in the far-ultraviolet (i.e. 
/l=1350-1750A) Afuv ^ 0.9, compatible with, but higher than 
the e xtinction est imated earl ier from SEP fits in the litera- 
ture (ISchaerer & PeU6ll2005l Ichary et alJboOSL Ipinlator et all 
20071k 

The Lya luminosity, = 1.9 x lO'*^ erg s"', is > 0.3% 
of the bolometric luminosity, which seems to foll ow the trend 



HCM 6A (z = 6.56) and GRASIL templates (Silva et al. IS 



noted for Lya blobs and S CUBA galaxies by iGeach et al 
(l2005b and lWebb et all (l2007h . 



3.3. Star formation rate 

The upper limit on the star formation rate (SFR) deriv ed from 
this FIR luminosity with the relation from Bell] (120031) assum- 
ing a Salpeter IMF from 0.1 to 100 Mq is SFR < 35 Mb yr" ' 



This u pper limit is a factor 8-2 lower than that of IWebb et al 



( 120071) and lCarilh et alj ( 120071) for Lya emitters at z ~ 5.7-6.6. 

Although the SFR depends strongly on the dust temper- 
ature (the SFR upper limit becomes 87M0yr"^ for rd=54K 
and 1 1 Mq yr ' for rd=18 K as indicated in Table|2]l this upper 
limit definitely excludes star formation rates observed in SMGs 
in excess of ~ 1 000 Mq yr" ' . Our upp er limit also exclud es the 
high SFR suggested for HCM 6 A bv lCharvetal] (l2005b . who 
attributed an apparent flux excess at 3.6 fim to He emission. 
The SFR limit is slightly lower but still compatible within the 
uncertainties with the estimates of the dust corrected SFR from 



Schaerer & Pellol (12005) listed also in Table|2] and based on the 
SED fits to UV-optical fluxes. Note that for comparison the lat- 
ter SFR values should be multiplied by a factor 2.55 to account 
for the different lower mass limit of the IMF (Miow - 1 M©). 

3.4. Star formation rate density 

The gravitational magnification allows us to probe a galaxy 
with a lower Lyo' luminosity, i.e. Li^ya,-2 x 10'*^ erg s"' instead 
of 7.5 X 10"^^ erg s ' fo r the average of the 2 galaxies discussed 



by I Webb et al. (12007 ). According to the luminosity function 
(LF) of iKashikawa et al. the cumulative number den- 

sity is at least 10 times higher at this luminosity and our upper 
limit on the SFR should therefore be more representative for 
the whole population of Lya emitters. 

To determine an upper limit for the infrared-measured SFR 
density (SFRD) we assume that the number density is given by 



the cumulative number density at LLyff=2 x 10"*^ erg s which 
is N(> L) ~ 1.x IQ-^ Mpc-^ For = 36 K, we then ob- 
tain SFRD(IR)< 2.4 x 10"^ Mq yr^' Mpc^^ with a possible 
uncertainty of a factor of ~ 3 due to the unkn own dust temper- 



ature (cf. Table |2]l. Integration over the LF of IKashikawa et al 
(12006.) down to 0.5L* (assuming a - -1.5, logL* - 42.6, and 
that Lijya has the same proportionality with SFR at all luminosi- 
ties) would imply an SFRD higher by a factor 2.2. 

Incidentally our up per limit for SFRD( IR) is very similar to 
the value estimated by IWebb et alj (l2007h . However, we note 



1 2 




1 2 3 4 5 

log(obs. wavelength) [/.im] 

Fig. 2. Comparison of the observed SED of HCM 6A (the data 
points are listed in Table[T]l wi th GRASIL templates for Arp 
220, M82, and NGC 6946 from lSilva et aP 119985 (from top to 
bottom in the IR, normalised arbitrarily at 4 ^m). The obser- 
vations in the near-IR and/or sub/millimeter range rule out Arp 
220 and M82 like templates, but are compatible with normal 
galaxies as NGC 6946. 

that these authors adopted a lower value for the cumulative 
number density drawn from the spectroscopically confirmed 
sample of Kashikawa et al. I (l2006h . whereas we use the LF from 
the complete z - 6.5 sample. The SFRD(IR) limit can be re- 
garded as a first estimate of the dust-obscured SFR density of 
Lya selected galaxies at z ~ 6.6. More observations are clearly 
needed to improve this estimate. 

For comparison, integration of the complete Lya LF, with 
the same assumptions as above, yields SFRD(Lya)=2.7 x lO"-' 
Mq yr"' Mpc"^ at z = 6.6 assuming a standard Lya SFR con- 
version of SFR(Lya)= 9. 1 X lO^^^'^LLya erg s ' . For well known 
reasons (dust and radiation transfer effects, partial IGM tran- 
mission) this should represent a lower limit of the SFRD of 
Lya emitters at this redshift. 

3.5. Comparison to empirical SEDs 
In Fig. 2 we compare the observed SED (cf. Table[T|) with semi- 



empirical templates from the GRASIL models of ISilva et al 
(1998). Shown are the templated of Arp220, a prototypical 
ULIRG, the prototypical starburst galaxy M82, and NGC 6946, 
an Scd galaxy. Figure 3 shows a comparison with the observed 
SEDs of low metallicity blue compact dwarf galaxies from 
Hunt et al. (2005). Clearly the observed millimeter to IR flux 
ratio excludes the ULIRG SED, and possibly also the dusty 
starburst (M82). This is in agreement with the moderate ex- 
tinction (Ay S 1) estimated from the fits to the (restframe) 
UV-optical observations. On the other hand, the SED of HCM 
6A is compatible with that of normal galaxies, like NGC 6946, 
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HCM 6A (z = 6.56) and BCDs (Hunt et al. 2005) 
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Fig. 3. Comparison of the observed SED of HCM 6 A with tem- 
plates of low me talHcity Blue Compact Dwarf galaxies from 
Hunt et al.l (l2005h normalised arbitrarily at 4 jum. The SED of 
BCDs, including the rest-frame UV part not shown here, are 
compatible with HCM 6A. 

and with the SEDs of blue compact dwarf galaxies (BCD). The 
present data does not allow us to distinguish between these two 
different types of SEDs, which mostly differ in their average 
dust temperature (dust being hotter in BCDs). However, HCM 
6A resembles more likely BCDs with higher surface brightness 
(or SFR surface density) than spiral galax ies. This is also found 
more generally for high-z Lyor emitters (IDow-Hygelund et al. 
2007h . In fact, local BCDs have SFR -0.2-3 Mq yr quite 
similar to the lower limit determined for HCM 6A. 

From the empirical SEDs shown in Figs. 1 and 2 we see that 
the millimeter flux may be close to a detection. In any case this 
object will be an interesting target for ALMA, which should 
be able to characterise the dust properties of this faint distant 
galaxy. 

4. Discussion and conclusions 

Quantifying the dust content in distant galaxies is important to 
constrain dust formation mechanisms in the early Universe and 
to determine cosmological properties such as the star formation 
rate density (SFRD) at high redshift. Strong gravitational lens- 
ing provides a unique access to the faintest observable sources 
probing the galaxy luminosity function down to more repre- 
sentative objects. Here we have used this effect to measure the 
millimeter emission of a faint lensed Lyo- emitter, suspected 
to suffer from a non-negligible amount of extinction in earlier 
studies. 

Our MAMBO-2 observations of the z = 6.56 galaxy HCM 
6A yield an upper limit of 1 .08 mJy (3cr) at 1 .2 mm. For typical 
dust parameters, this translates into upper limits of ^ 5 x 1 0^ M© 
for the dust mass, ^ 2 x 10" for the infrared luminosity, and 



S 35 Mq yr ' for the IR star formation rate (SFR). Thanks to 
strong gravitational lensing and to excellent weather conditions 
these limits improve recent sub-mm measurements of two z ~ 
6.5 Lya emitters by approximately one order of magnitude. 

Assu ming a dust production of ~ O.I Mq per type II su- 
pernova (IBianchi & Schneideii 120071) a maximum age of ~ 
1 /2f// =s 4 X 1 0^ yr, would require a constant SFR ~ 60 Mq yr" ' 
(for Miow - 1 Mq) to produce 5x10^ Mq of dust. Hence our 
results are marginally consistent with current estimates for the 
dust production. Reducing the large age uncertainty of HCM 
6A, and lowering the upper limit for its dust mass will allow to 
place more severe constraints on dust formation scenarios. 

Our non-detection of dust does not exclude the presence of 
dust in HCM 6A. In fact the SFR(IR) upper limit of S 35 Mq 
yr"', with a typical uncertainty of a factor 3 in either direc- 
tion, comes close to but remains co mpatible with the UV SFH 
of ~ 28-105 Mo yr obtained by iSchaerer & PeUd (l2005 in 
for an extinction of Ay ~ 1. Its is tempting then to speculate 
that the IR emission from this object could soon be detected 
with somewhat deeper sub/millimeter observations. However, 
another explanation for the apparently red SED of the Lya 
emitter HCM6A may also be a composite stellar population 
with little or no dust, as discussed by lSchaerer & Pellol (l2005l) . 
Deeper sub/millimeter observations of lensed high-z galaxies 
including HCM6A should therefore soon be able to match the 
constraints from shorter wavelength observations and provide 
more detailed constraints on dust in the early Universe. 

Thanks to gravitational lensing the Lyo' emitter observed 
here coiTesponds, with LLyir = 1.9 x 10"*^ erg s"', to a ~ 0.5L* 
ga laxy according to the z ~ 6.5 Lyo- luminosity function (LF) 
of Kashikawa et al. I (l2006h . at the faintest limits of their blank 
field survey limits. Using our millimeter non-detection and 
their LF we determine an upper limit for the infrared-measured 
SFRD of < 2.4 X 10"^ Mq yr"' Mpc"^ with an uncertainty of a 
factor 3 in either direction. This can be regarded as an estimate 
of the dust-obscured SFR density of Lya selected galaxies at 
z ~ 6.6. More observations will be needed to improve this esti- 
mate. Gravitational lensing and new upcoming facilities should 
allow significant progress in the near future. 
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